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I. INTRODUCTION 

One of the main interests in hadronic physics is to con- 
struct effective, low energy approximations to QCD and 
to find methods enabling to treat the theory in its non- 
perturbative domain. Over the years a number of non- 
relativistic or semi-relativistic models describing quarks 
in hadronic bound states have been proposed, however, 
the full effect of quark-antiquark pairs has never been 
considered. This is because such effects lead to a many- 
body problem that can only be treated in some approxi- 
mation. In hadronic models such approximations are typ- 
ically driven by phenomenological considerations rather 
then QCD itself. Here, instead we construct a simple ef- 
fective quark Hamiltonian guided by QCD and look for 
various classes of solutions which may represent hadronic 
states. 

A possible connection between QCD and an effec- 
tive quark model Hamiltonian was presented in Ref. 
A confinement scenario and dynamical chiral symmetry 
breaking were discussed, however, even in the chiraly bro- 
ken phase the effective Hamiltonian remains complicated, 
i.e it contains interactions among an infinite number of 
particles. In Q following certain assumptions, a simple 
solution for the SU{2) case was presented. In particu- 
lar quarks were assumed to be confined in single s-wave 
orbitals in a finite volume and no dynamical contribu- 
tions from gluons was considered. The resulting effec- 
tive Hamiltonian allows for analytic solutions with low 
energies saturated by color-neutral physical states and 
colored states shifted to arbitrarily high energies. That 
study was followed by [| were the limitation on the num- 
ber of quark orbitals was relaxed and more analytic so- 



lutions were found. 

In this work we extend the method of Q to include 
all orbital levels in the quark sector. The color part is 
now extended to SU{3), while in the flavor sector we 
keep the approximately mass-degenerate u and d quarks. 
As in dynamical contributions of gluons are not con- 
sidered, and the system is confined to a finite volume. 
Under these approximations we show how QCD can be 
solved nearly analytically, if the potential interaction of 
the QCD Hamiltonian is replaced by a spatially averaged 
constant interaction. 

The paper is organized as follows: in Section II the 
model space and the Hamiltonian are defined, In Section 
HI we study the analytic solutions for two- and three- 
orbital level systems, in the chiral limit. This is followed 
by the consideration of an arbitrary number of levels, 
and the quark mass term. We discuss why in that case it 
is impossible to find simple analytic solutions, unless an 
orthogonal transformation and a BCS-type of mean field 
scheme is adopted. In section IV we discuss the possible 
emergence of meson and baryon spectra that result from 
the lowest s- and p-quark orbitals. Conclusions arc drawn 
in Section V. 



II. THE HAMILTONIAN AND ITS MODEL 
SPACE 

The construction of an effective, QCD-inspired Hamil- 
tonian that approximates its low energy spectrum is 
based on the following assumptions. Due to confinement 
the domain of fields is expected to be restricted to a finite 
volume in space, where individual hadrons are located. 
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FIG. 1: Illustration of the effect of the Hamiltonian on the 
basis states. The operators, which appear in the kinetic en- 
ergy, can be divided into columns of operators, related to a 
given spin j, which commute with those of different columns. 
Only the positive energy states are plotted. The ellipses indi- 
cate which states are connected through the interaction terms 
appearing in the kinetic energy. Columns of different j are 
separated by a vertical dashed line. 



This is achieved by using the 3-dimcnsional harmonic 
oscillator for the orbital states. The confinement scale 
is then related to the width of the harmonic oscillator 
wave functions, 7 De-confinement transition can then be 
studied as a limit when 7^0 since in this basis for a 
quark to probe large distances would require mixing with 
a large number of excited states. Alternatively one can 
use a 3-dimensional box, but, as we shall show later on, 
the oscillator representation is mathematically easier to 
handle. Gluons arc not taken into account dynamically, 
rather they are accontcd for by a static potential and 
the effects of a confining potential arise from the use of 
a confining harmonic oscillator. Finally the potential is 
smeared over the hadron scale and effectively replaced by 
a constant. 

Most of the restrictions imposed in 0, Q are now re- 
moved. Therein we restricted quarks to be in SU (2) color 
but, as we shall see here, this restriction can be easily 
eliminated. We shall also include all the orbital levels, 
by following the method developed in 0]. 

In Fig. [T]the model space for the quarks is schemati- 
cally presented. Only the orbital states at positive energy 
are depicted, while vertical dashed lines separate columns 
of equal total spin j. As it will turn out, the Hamilto- 
nian can be written in terms of operators which act only 
within each column. This will greatly simplify the task 
of finding an analytic procedure. 

The fcrmion creation operator for a quark in the level 
a (a = ±i denotes the upper/lower Dirac level), with 
angular momentum I coupled with quark spin- 1/2 to to- 
tal spin j, is given by 



ma 

The remaining indices: N, A, c and / refer to the oscil- 
lator number, spin, color and flavor index, respectively. 
The (Zm, icT | jA) is an SU{2) Clebsch-Gordan coeffi- 
cient. The quarks are in the spin and flavor-spin repre- 
sentation i , while for the color part they are in the triplet 
irreducible representation. The creation operator on the 
right hand side corresponds to a spin-orbit decoupled ba- 
sis. The operator with upper-level indices is defined in 
the standard way, (see Eq. [T]), 

^-a(N,l.l)j-X.c-f (2) 

Here and in the following the capital index C denotes 
globally the color part. F gives the flavor-spin. The 
phase convention is taken from 0] (with a corrected sign 
in front of T^). The lower-case c denotes the magnetic 
quantum numbers of the color part and is a short hand 
notation for (Yc, Tc, Tcz), with Yc being the color hy- 
percharge, Tq the color-isospin and Tcz its third com- 
ponent. The bar over the index c refers to change from 
c = {Yc,Tc,Tcz) to c = (,-Yc,Tc,-Tcz) (the conju- 
gate index). In the following we abbreviate (— l)~f'"'"^'^ 
as (-1)^^=. 

When the index a has the value a = +^ (upper level), 
the operators b are replaced by a, and are referred to as 
quark operators. In contrast, when a = — ^ we define 
them as d and as the anti-quark operators. The b cre- 
ation/annihilation operators are then converted into the 
d annihilation/creation operators. This notation exploits 
the Dirac formulation for fermions, where a fermion in an 
upper level is denoted as a particle while a fermion hole 
in the lower level is denoted as an anti-particle. These 
and their conjugate will be used to express the effec- 
tive Hamiltonian. The kinetic, Dirac energy is originally 
given by 

K = J dxijj\x)[-iV ■ a + (3mo]i>{x). (3) 

A mass term has been included, in contrast with Q and 
the potential term is given by 

V = J dxdyil)\x)T''ip{x)V{\x-y\)ip''{y)T''xp{y). 

(4) 

The static potential V{\x — y\) simulates the gluon-quark 
interaction. The origin of this expression is the Faddeev- 
Popov term of the QCD-Hamiltonian The latter is 
essentially a color-color, Casimir interaction and when 
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averaged over the confined hadronic bag it is replaced 
by a constant Vq, and becomes exactly the total color 
operator This can easily be seen by noting that 

/ dxil>\x)T'^rp{x) represents the color charge and a con- 
stant potential separates the integrals over x and y. 

In what follows, we will first consider the potential in- 
teraction, because it will result in a very simple spectrum 
which can later be used to diagonalize the hopping inter- 
actions brought by the kinetic energy. 

To do so, we expand the fcrmion fields i/j^ and xp in 
this, complete particle basis in terms of the decoupled 
fermion creation and annihilation operators, 

xpix) = ( 

with 

Nlm 

■^2(3;>Cr:C,/) = '^b'^'^"''"'"'^ RNL{\x\)Yirn{n^)x<T- 
Nlm 

(5) 

Here fix denotes the angular components 9^ and (j)^- of 
the postion vector x. The fields have components in spin, 
flavor and color. In Eqs. (l3|),(|4]) for the kinetic and po- 
tential energy, sum over the spin, color and flavor indices 
is implicit. 

A. The potential 

Substituting into Eq. ^ the fcrmion fields, using 
Eq. ([5|) and transforming to the coupled fcrmion creation 



and annihilation operators as defined in Eq. (H)) leads to 
the following form for the potential interaction (see Ap- 
pendix A for details) 



V^Y. v/(2ji + l)(2j3 + l)m , h , k , L) 

X 

(jiAi, j2A2|LA/)((l,0)ci, (0, 1)C2|(1, l)c)i] 

a',{N3,l3,^)j3,\3,C3f' (7V4,i4,5)i4A4,C4 

(J3A3, J4A4IL - M)((1,0)C3, (0, 1)C4|(1, l)c)i] 

(6) 



where X stands for the collection of: 
Ni, ji, Xi, li, Ci, a, /, a', /', L, M, c quantum numbers. 
The SU{3) Clebsch-Gordan coefficients carry a multi- 
plicity sub- index, 1, because the fermion creation and 
annihilation operators are coupled to the generators 
of the color-S'?7(3). The coupling corresponds to the 
(1, 0) (g) (0, 1) = (0, 0) + (1, 1) decomposition of the prod- 
uct of two quark representations in which no multiplicity 
appears. The convention is taken from 0, @, 0] (with 
a corrected sign in the 5'[/(3)-phase) . As before, the 
indices a, a' are used to label pseudo-spin components, 
and Zi, /, Ci, ji, Xi label orbital angular momentum, fiavor 
and color components, total spin and its projection, 
respectively. The quantities V{Ni,li,ji, L)^s define the 
intensity of each component of the interaction in Eq. 31 
and, (before space-averaging) are given by 



V{N,,k,j.,L) - ^^^^^^^ (-1)^'- + ^+^"^+^^ I |a;|2d|a;||y|2d|y|i?jv,;,(|a;|)i?^2;2(|a;|)i?iV3;3(|y|)^A^4;4(|y|) 
X /\(cos.)P,(cos.)y(|.|, |y|,cos.)ntr Vi^^k + 1)(2,. + mMOmUMOm / /, ^3 i 



V(2ji + l)(2j3 + l) 



(7) 



The factor | originates in the triple-reduced matrix ele- 
ments of the color-operator, as explained in Appendix A, 
and 9 is the polar angle angle between vectors x and y. 



For a space averaged V(|x— y| Vq potential, the only 
non vanishing contribution comes from L = 0, which as 
a consequence separates the integral over x and y and 
leads to V{NiJi,ji,L) to 



Vq 

V{Ni,li,ji,L = 0) = -^5L,Q5NiN2^jij2^iii2^N3Ni5j^j^5i^i^ 

(8) 

Thus for the potential term in the Hamiltonian we fi- 
nally have. 



V = ^C2(5[/(3)), 



(9) 
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where C2{SU{3) is the second order Casimh operator of 
S'f7(3)-color, given by 



N 



C1C2 



ci 

C2 



-N 



"a(Ar;,i)jAci/" 

aNljXf 



,a{Nl,^)jXcf 



aNljXcf 



(10) 



The C^_^ are the generators of the SU{2>)-co\ot group. 
Its eigenvalue is given by Q , 



(11) 



with (Ac,/^c) defining the S'C/(3) irrcp. Thus, the po- 
tential part simply separates colored from non-colored 
states. The second-order Casimir operator is equivalent 
to the color-spin in S'{/(2)-color discussed in 0. 



B. The kinetic energy 

We first consider the mass term in Eq. ([3]). It leads 
to an operator proportional to the sum of the quark and 
antiquark number operators: 



mo(nq + Uq) + mo Eat/jAc/ 1 



-mo 



_ht l,-^(JVj + i,i)jAc/\ 

"-iCATj + i.ib-Ac/" ) 

l^N'jXcf ^Oi(jV'j-i,ibAc/'' 

''-i(iv'j-i,i)jAc/'' ; ■ 



(12) 



The terms with a = 5 count the number of quarks in the 
upper level, while the ones with a = — i together with 
the last term count the number of holes, thus the number 
of anti-quarks, in the lower level. The last term in the 
first line is a constant and may be skipped, if convenient. 

The rig = b\ b^" = ala" and rig = b'^'^bl, = d^'^d^ 

■2^ 2^ 
are the quark and antiquark number operators, respec- 
tively. The index is a short-hand notation for all the 
indices which label individual creation and annihilation 
operators. 

The momentum dependent part of the kinetic energy, 
when expressed it in terms of the fermion creation and 
annihilation operators, is given by (see Appendix for 
derivation) , 



K 



where. 



(13) 



N^j + 



3\ 2 
2 



"i(JVj + i,ibAc/ 



^-i(Af+lj-i,i)iAc/ 



^-i(JV-l,i-i,i)iAc/ ^ 



N-J + 



3\ 2 
2 



"i(W+lj-i,i)jAc/" 



^- - v^E E E 

j N=] + ^ ^cf 



-i(Afj + i,i)jAc/ 



b' . 



(Af-lj-i,i)iAc/ 



^i(Ar-lj-i,i)jAc/ _^ 
^i(Wj + i,i)iAc/ 



^i(W+lj-i,ibAc/ 



N-J + 



3\ 2 
2 



/,t .i(iVj + i,i)iAc/ 

-i(7V+lj-i,i)jAc/" 



(14) 



The X+ operator shifts quarks from the lower {a — —1/2) to the upper (a — +1/2) level, while K- does 
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the opposite. The tilde is used to distinguish the above 
operators from actual generators of a SU (2) algebra. 

The mass term will be skipped in most of the cases dis- 
cussed here for the following reason: the kinetic energy 
part without the mass term can be identified as being 
proportional to a component of a generator on an SU (2) 
algebra. Thus, it can be diagonalized exactly. The mass 
term, however, does not commute with the K^^ opera- 
tors, which determine the momentum dependence of ki- 
netic energy and destroys exact diagonalizability. That 
they do not commute can also be seen noting that, for 
example, the operator creates a particle-hole pair 
which is equivalent, as just noted, to raise the number of 
quarks and antiquarks. Nevertheless, when we treat the 
complete problem, though the SU{2) structure is lost, 
we can still diagonalize the kinetic energy part, using the 
BCS formalism M. 



III. ANALYTIC AND SEMI-ANALYTIC 
SOLUTIONS 

In what follows we shall consider two, three, and fi- 
nally an arbitrary number of single-quark orbital levels. 
The reason for choosing this path is to investigate the 
analytical properties of the solutions and if any features 
in the small basis approximation can be generalized to 
the case when any number of quark levels is allowed. In 
the two- and three-level case, the mass term is neglected, 
otherwise no simple solution can be obtained. The mass 
term will be included when all orbital levels are taken 
into account, since as we show below in this limit the 
mass term can be treated together with the kinetic en- 
ergy within the BCS formalism. 



A. The two level-systems 

This refers to case when, for a given j two orbital levels 
are considered, one with N oscillator quanta and with 
orbital angular momentum I = j + ^ and the another 
with either N' = {N — 1) or (N + 1) quanta and orbital 
angular momentum I = j — ^. For example, for = 1 
and j = \ this corresponds to the lowest p- plus the 
lowest or next-to-lowest s-orbital. 

The momentum dependent part of the kinetic energy 
can be related to K± operators. These operators are 
proportional to Km, when restricted to two orbital levels, 
and are of the form 



000 

000 
000 



000 



(15) 



The spin-color-flavor coupling denoted by " (g)" is defined 

as 



(16) 



Here and T stand for the quantum numbers dent- 
ing combined representation of the spin, color and flavor, 
and /Jk, M denote shortly the magnetic quantum num- 
bers, e.g., fj, = Xcf. The expression in Eq. (fT5|) can be 
easily derived from Eq. (|14p by noting that the coupling 
to spin-color-flavor singlet implies a contraction of all in- 
dices. The factor y/2 ■ 8 ■ (2j + I), is the square root of 
the multiplicities of the flavor, color and spin represen- 
tation, respectively, and yields the correct normalization 
of the singlet representation. 

The first index in the creation and annihilation op- 
erators refers to the pseudo-spin component (i) for the 
upper level, and (—5) for the lower level. Thus when 
1 = 1, the kinetic energy couples quarks in the s- level 
with quarks in the p-level, of total spin j = \- When 
I = 2, quarks in the p-level, with total spin J = §, are 
coupled by the kinetic term to quarks in the d-level, also 
in the total spin |. The orbital values I are given by 
j ± i. The radial number, TV, starts from j + ^ and ac- 
quires odd (even) values, while N' has even (odd) values 
only and satisfies N' = N — 1 or N + 1. This selection 
rule is strictly obeyed only for the harmonic oscillator. 
Note, that the Km operators are coupled, for each given 
spin j, to total color zero. Thus, they commute with the 
total color-spin operator of each column j. 

The if-operators satisfy the commutation relations 



^0 i^i 



±K 



jN 



(17) 



i.e., for each combination of N and j they form an 
SU^^{2) pseudo-spin group. Furthermore operators with 



different spins commute. The operator Kq^ is given by 



K 



jN 



K 



■jN 



/2~8v/27TT 



{N',l,^)j'^^^{N,l+h^)j 



000 

000 
000 

000 



= V2^^/2j + 1 



- In- 



i.e. as a half of the difference between the number of 
quarks in the upper and lower level. The quarks are cou- 
pled to total spin j. Note, that the pseudo-spin group 
SU^^ {2) defined by the Km operators is different from 
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the pseudo-spin group in the original definition, which 
has been given with reference to lowering and raising op- 
erators within the same orbital. The K,n operators, de- 
fined here, raise and lower quarks from one orbital to a 
different one. 

The relation to the operators in Eq. p4|) is obtained, 
by either (a) restricting to the combination of the orbital 
Nj=j + l with the orbital {N - = j - ^, or (b) 
to the combination of the orbital N,l = J + 5 and the 
orbital (N + 1),/ = j — ^. The relation of the kinetic 
energy term and the generators of the SU{2) algebra are 
is then given by 

for the case (a) 

= /^^^^^/^r-Wi^ (19) 
and for the case (b) 

^ ^EiMk;" , A^^^Ki". (20) 

The operators commute with the total color oper- 
ator in each column j and, thus, also with the total color 
operator. This is because for each given j the they arc 
coupled to color zero. As a consequence, the Hamiltonian 
can be diagonalized analytically. It is given by 

H'" = A.N, [k^ + K=^) + ^C2(5C/(3)),(21) 

with A.Mj (K=a,b), defined in Eqs. (HH) and (HO]). The 
expression in the parenthesis of the first term is just K-fJ^ , 
whose eigenvalues are known. The kinetic part is just 
then 2Ka;. Thus the eigenvalues of this Hamiltonian are 



+y (Ac + AcMc + yuc + 3Ac + i^ic), (22) 

where A'lj is the projection of the pseudo-spin operator, 
within the pseudo-spin J, onto the x-axis and (AcMc) 
define the irreducible representations (irrep) the SU (3) 
color group [1, 01- When color is zero, the irrep is given 
by (0, 0), while for a color octet it is (1, 1), etc. We used 
the definition of the second order Casimir operator as 
given in 0, Q . For color less states this reduces to the 
simple formula 2Af^pfjMj. 

Fig.[l]illustrates the effect of the kinetic term discussed 
above, without showing the negative energy levels The 
negative energy levels are just a copy of the levels at pos- 
itive energy inverted to negative energy of what is shown 
in Fig. [1] Each of the ellipses represents a combination 
of orbitals for which the analytic solution given above in 
Eq. applies. The orbitals, which the ellipses connect, 
are not only the ones at positive energy, as shown in the 
figure, but also represent the connection to the negative 
energy levels. 

B. The three-level system 

The three-level system consists of the following levels: 
We add to each level, given by j, N,l — j + ^, two levels 
with I = j—^ and (A'^ — 1) and (A^ -|- 1) oscillator quanta. 
The kinetic energy for this three-level system becomes. 



k = k\ +k_ 

where 



K 



jN 



K 



Ac/ 



A^-J + 



l(Af+lj-i,i)jAc/ 



1^ 2 ) ''^(JV-lj-i, lb-Ac/'' 



Ac/ 



"-i(Wj + i,i)jAc/" 



i(Ar+lj-i,i)jAc/ 



/,t /,-^(JVj + i,ibAc/ 



''-i(7Vj + i4)jAc/'' + 



"-i(7V-Hj-i,i)iAc/ 
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-i(7V-l,i-l,l)jAc/ 



(24) 



Note, that this is exactly the same expression as in 
Eq. (|14[) but without summations over j and N. 
Defining 



K 



jN _ 



jN 



K 



T]K 

„~3N 



(25) 



it can be verified that, by choosing the factor rj appro- 
priately, the operators K^j^ , just like in the two-level 
case, show the standard form of a SU{2) algebra, such 



that 



K 



= Ki^ and 



K 



±K 



jN 



(see Appendix D for details). In the three- level case the 
result then resembles the one in Eq. (|2ip including ex- 
pression for the energy given by Eq. (|22p. This comes to 
a complete surprise because it is not obvious at all that 
an SU(2) structure is contained in the relativistic kinetic 
term in this quasi-particlc basis. 

The fact that we found even in such a complicated 
system an analytic solution, is an evidence that probably 
a complete analytical treatment, of the whole problem, 
i.e., including all orbital levels, might be possible, save 
a simple numerical solution of a set of equations. We 
discuss this case in the following section. 



C. An arbitrary number of levels 

The kinetic energy that is not restricted to operate 
in a specific subspace of single particle orbitals, is given 
by Eq. p4p . The two parts of this kinetic energy, the 
one which moves quarks to the upper level and the one 
which moves quarks to the lower level, do not satisfy a 
SU{2) algebra any more. We will also add a mass term, 
which by itself destroys any SU (2) structure once present 
(see discussion in II. B). Nevertheless, we can exploit the 
structure encountered in the two- and three- level case to 
simplify the problem. 

Noting that the total kinetic energy is given by 



K 



(26) 



we can select one particular value j and solve the problem 
for it. The complete solution is the sum of the solutions 
for all j. We introduce a cut-off value (n) for N, which 
can be chosen arbitrarily. Therefore, 



n.AN=2 min(n.Ar+l) 

N=j+^ N'=N-1 Xcf 
n,AN=2 max(n,Af+l) 

N=j + i N'=N-1 Xcf 

n 

+"^0 2^ 2^ ^»i(jVj + i4)jAc/'' 0^1(jVj + i,i)jAc/'' ) 

N=j+^ Xcf 

n 

+mn V Vffot bK^^'j-i^b'^c/ _ .t ,_i(Ar'j-i,i)jAc/\ ,27) 

+™o 2^ 2^[''^iN',j-^Si^-f ^-i(w'j-i4)jAc/'' )■ ^^^) 

Ar'=j-i Xcf 

I 



In the last two terms, the step size of N and N' is 2. We 
have skipped a trivial constant at the end of Eq. ([27|) . 
which will be included in the BCS formalism that we 



treat later. The factors k-j^j^, can be read from Eq. ([14]). 
The kinetic energy can now be rewritten as 



AN=2,n min(n,Ar+l) 

1^0-,n) _ V V ViP \ fh^ /,-i(JV'j-i,ibAc/ ,t hi(N',j-i,i)jXcf\ 

N=j+^ N'=N-1 Xcf 
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AN=2.n min(n,Af+l) 

Ar=j + i JV'=JV-1 Ac/ 
AAr=2,n 



Ar=j + i Ac/ 

"^0 2^ 2^ ^''i(W'j-i,ib-Ac/'' ''-iC 

7V' = ,-i Ac/ 



^5-5(JV'j-i^bAc/^ 



(28) 



while the mass term is unaffected. Next, we apply a 
unitary transformation 



±l{N,j + i,i)jXcf 



-*iArfc"±i(fe.j + i,i)jAc/ 



^±i(JV'j-i,i)jAc/ - I]^jV9^±i 



The new operators have to satisfy also the fermion anti- 
commutation rules, which imposes a constriction on the 
a's and /3's. In particular for the operators with orbital 
spin I = j + ^ we have 



/j,t J,±5(^2j + i,i)jAc2/2\ 

\"±i(WlJ + i,i)jAci/2'" / 

= SNiN2^CiC2Sf\f2 

=^ J2k(^'jNiki^jN2k2 = SNiN2^kik2 ■ (30) 

while for I = j — ^, 

\°±i(iv;j-i,ib-Ac/'" / 



— SN[N^SciC2Sfif2 



(31) 



A few words are necessary, concerning the variables 
ajATfc and f3jN>q- We cut the space at a given N = 
-^max = for J + ^, such that the maximal N' equals 
Nniax ~ 1- Choosing a given integer n as the maximal 
number for A^, the range of fc and q is from 1 to Uj (see 
Fig. [2]), where is the number of levels for ^ = j ± i. 
For example, in the two level system and j = i, we con- 
sider one p orbital (orbital spin j + 5 and = 1) and 



one s orbital (orbital spin j 



\ and A' 



0). In this 



case, Uj — 1 and k = q = 1. Using the symmetry of the 
matrices, we get for ofj^vfe a total of nj{nj + l)/2 linear 
independent elements, and also for Pjj^g- In order to de- 
termine all of them, we need njijij + 1) conditions, as we 
shall explain later on. Applying these transformations, 
the kinetic energy reads 



T^3 X^X^ZO fu'< T-i{qJ-hi)3^'^f , rt ^i(9j-i,i)jAc/\ 

^ = Z^2^''kq[b^ik.j+^,^)jxcfb +b_i(fe_^.+ i_i)^-A,;6 1 

kq Xcf ^ ^ 

, Y^Y^7j r-i(fej+|,^)jAc/ ^i(fcj+i4)jAc/\ 

■^2^1^^kq\b\{q,,-\,\)jXcsb +b_l(g,,_l^l)jAc/b I 

kq Xcf ^ ^ 

, V^V^ /rt ^:^(fcj + i,i)jAc/ „i(fej+i,i)jAc/\ 

+ 2^2^"'0,fcj + i I + -''-i(fej + i,ib-Ac/b j 

fc Ac/ ^ ^ 

, Y-Y^ /rt ^i(g,j-i,i)jAc/ ^-i(g,j-i,i)jAc/\ 

9 Ac/ ^ J 

I 

with 

Tl-l 
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variables 



variables 



n-1 



FIG. 2: Structure of the matrices a j (3. The rij gives the 
dimension of the matrix. Though A'^ goes from j + \ up to n, 
the maximal number of A'^, considering that A'^ makes steps 
in two, the number of rows are exactly rij. the relation of rij 
to n is given by = i (n — (j + i)) + 1. 

n 

mo,fcj+i = ^ mo\ajNk?- (33) 

A'=l,3,5,.. 

I 



and the factors fc^^ are given by 

AN=2,n min(n,A''+l) 

= E E \^NN'W)NkP,N',. (34) 

Ar=j + i N'=N-1 

In order that the kinetic energy, for a given j value, is 
diagonal in the orbital index, we require 

^ for fc 7^ g. (35) 

These leads to nj(nj — 1) conditions which together 
with the 2nj normalization conditions, | ajNk |^ = 
1 and I PjN'q ^ — 1, leads to nj{nj + 1) equations 
required to determine all variables UjNk and PjN'q- 

It shows that at the end, when all restrictions are ful- 
filled, we arrive at new orbitals labeled by the index 
fc = 1, 2, .., rij, which are divided into particles and anti- 
particles. Note that, when only the diagonal components 
of k^.^ are different from zero, the kinetic energy con- 
tributes to Tij new orbitals for each orbital angular mo- 
mentum I = j + and / — j ~ \- -By introducing the 
index k, the final form of the kinetic energy reads 



k Xcf 



EE^ife 

k Xcf 



j-ll)jXcfb 



■b_i^k,j-i^}jxcfb 



k Xcf ^ 
k Xcf ^ 



b — o_ 1 



.i(fc,i-i,i)jAc/ 



b-i(k,j-^,^)jXcfb 



-l{k,j-^,^)jXcf 



(36) 



The number-operators of quarks and anti-quarks (uq 
and Hq) are invariant under the unitary transformation 
described above. They just transform into the new oper- 
ators with that also represent number operators of quarks g2('=j'+2>2)j-^'=/ = C- 1 Jj^'^'i^^'^f ^ g . 1 ^Hk,j}Xcf 



or anti-quarks, the latter represented by holes in the 



L ^ w . I 1 ^ ± i.Uj . 



lower levels. The same is true for the mass term - it l lik,j+l,l)j^cf ^ •{k,j)Xcf i'\(k,j)Xcf 



does not change under this transformation. 



Up to this point, however, fermions with the orbital ''-i(fej+i,i)jAc/ = *j+i,fc^]- i (fcj)Ac/ + S+|:fc'^i+i(fcj)Ac/ 

angular momentum / j ± i are still mixed-in by the i , 1)jAc/ _ Ak,j)Xcf .t(fcj)Ac/ 

kinetic energy. The final diagonalization will be achieved, " ~ '^j+^,k"j-i ^ 
applying a BCS transformation Q. The transformation ^i(^k.j-^.^)jXcf {k j)Xcf ](k j)Xcf 

is given by b' =Sj_^j.b^^^ +Cj_i f,d._[ 
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Oi(fej-i,i)jAc/ 
Oi(fej-i,i)jAc/ 



ubov angle. Applying this rotation to the kinetic energy 



^ term and including the mass term, leads to (Irom now 

Sj-i,fc^]+i(fc.j)Ac/ +S-^,fe'^j-^(fej)Ac/-On we skip the upper index j in fc^^, because from the 
^ , , context it is clear that we work in a fixed column denoted 



with c = cos(0) and s = sin(6'). and 6 being the Bogol- 



^^BCS 



J dx'ip\x){—ioL ■ V + f3mo)'i{){x) 



Ac/ fc 



6^ 1,, .,, ,6. 



,(kj)\cf 



aj-i(fei)Ac/ 

j,t ;,(fej)Ac/ 
j + i(fcj)Ac/"j + i 

^t(fcj)Ac/ , 



- (2 •2-3) 



3 ^ ^ 1 



L,fc) + ™0,fcj + i(c'+i,fe-S-_i,fc) 



+ 2kkk i,fcCj-_i,fe + Sj+i,fcCj+i^fe 



. J-^(fej)Ac/^-+i + S+i(fcj)Ac/'^j_ 1 
Lt jt(fci)Ac/ , , (fcj)Ac/l 1 



(fej)Ac/ 



r 



(38) 



The kinetic energy operator is diagonalized, requiring operators vanish, which leads to the gap equation that 
that the terms quadratic in the annihilation and creation determines the Bogolubov angle, 

I 



kkkic^i 



= 
= 0. 



r 



(39) 



The diagonalized new operator, for a given column j, is finally given by, 
I 



"■BCS 



/d.^t(,)(_,,.V + /3^„)^(. 

Ac/ k 



2 2 



^-i(fej)Ac/"j-i 



Jfej)Ac/ 



2 2 



jt(fcj)Ac/^ 



3j-i(fci)Ac/ 
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+ 



^t(fcj)Ac/^ 
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I 



(40) 



Introducing a notation AS to represent the 
Sac/ ^XcfB^'^^ it can be rewritten as. 

A: 

(41) 

This aUows us to determine the spectrum of such a 
model including excited states, which are to be associated 
with hadrons. In particular meson states are constructed 
from 



](0,0}cF 

J(0,0)c-F 
jiO,0)cF 

j{0,0)cF 



(42) 



where the j refers to the total spin, (0, 0)c indicates the 
overall color singlet and F gives the total meson flavor. 
The quantity /i is a short-hand notation for all mag- 
netic quantum numbers. The colored states are shifted 
to larger energies due to the second order color operator. 
The structure of the meson spectrum is already quite in- 
volved. Examining (|42|) and apply to these states the 
energy operator of ()4ip . we obtain the meson spectrum. 
Because in (|4ip no explicit dependence on the coupling 
to spin and flavor appears but only on the new orbital 
indices k and k' , states with the same spin and flavor 
are degenerate in energy, fixing k and k'. Only when 
the orbital indices k and/or k' change, also the energy 
changes. For example, we can expect that the lowest 
states {k = k' = 1) with flavor-spin 1 and multiplets 
with j ~ and 1 are degenerate. 

For the baryon spectrum, the states which correspond 
to three quark operators (b^ or d)) should be coupled 
to the desired flavor and spin. The nucleon is then de- 
scribed by the three quark operators which correspond 



to the lowest energy. Also here, P ~ j = mul- 
tiplet, containing the nucleons, will be degenerate with 
the ^ = f J J = I multiplet, containing the A mesons. 

Each quark operator is a complicated function of the 
bare quark operators of Eq. ([5]). Thus, in therms of the 
bare operators hadrons states contain a sea of quark- 
antiquark pairs in the background. 

As it can be seen, the kinetic energy operator, includ- 
ing the mass term, can be diagonalized exactly. In the 
next subsection we apply the procedure to a four level 
system and show how the bilinear equations, involving 
the a and /3 matrix elements, can be solved. 



IV. A TEST CASE: THE TREATMENT OF 
FOUR LEVELS 

We consider the four level case, with j ^ ^, including 
the mass term. We take the first two p {N = 1,3) and 
first two s {N' = 0, 2) orbital levels. The coefficients k^.^ 
and ajNk, PjN'q satisfy the relations 



N'q\ 



= 1 



E 



ajNk] 



(43) 
(44) 

(45) 



where 



^kq 



min(3,7V+l) 

E E KN'\»*Nkf3N'q (46) 



Af=1.3 N'=N-1 



and the factors | kj^j^, \ are given by 




N-j + 



0N',N + 1 



7i 0N'.N-1 



(47) 



This expression shows that = 1 connects to N' = 
0, 2 and A^ = 3 connects to A^' = 2, which is due to the 
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Solution 


Oil 


ai2 


Q32 




/302 


/322 


kii{GeV) 


fc22(GeV) 


1 


-0.53452 


0.84515 


-0.53452 





1 





0.18092 


0.23357 


2 


0.57930 


-0.81511 


0.57930 


-0.10050 


-0.99493 


-0.10050 


0.14320 


0.22184 


3 


-0.99004 


-0.14072 


0.99004 


-0.99276 


0.12010 


-0.99276 


0.30450 


0.32383 



TABLE I; Relevant solutions for q's, /3's, for the the system of four levels. The total spin of the level is j = ^, i.e. the first 
two s and p orbitals are considered. 



use of the harmonic oscillator basis. The scale ^/j for 
this calculation is taken to be of the order of l/m~^. 
From here on we drop the upper index j in fc^jv'- "^^^^ 
conditions that need to be satisfied are, 

"11 + "12 = 1 
"12 +"§2 = 1 
/?o'i+/3o'2 = l 



|fcio|aii/3o2 + |A;i2|aii/322 + |A:32|ai2/322 = 
|fcio|ai2/3oi + |A:i2|ai2/3o2 + \k32\a32i302 = 0, (48) 

where a's and /3's are real. These equations are solved 
according to The numerical solution of these equa- 
tions shows the existence of several sets of solutions (an, 
"12, "32, /3oi, /3o2, ^22)- These solutions produce re- 
peated values for the coefficients kkk- Finally, there are 
three sets of solutions that produce three different k^k- 
These sets are described in Table I. 



In the transformed basis the kinetic energy including 
the mass term can be written as 



K%ug = X] X] { \^kkksp^kcp^k + rno^k,sc%^k ~ "io,fc, psp j. 



'S{k3)\cf''S 



Xcf k 

^2kkkSs,kCs,k + mo^k,scl^k - ™o,fc,pSs,fc 
'2'kkkSs.kCs,k + rno^k,pc's,k ~ rno^k,ssl,k 
2kkkSp^kCp,k + "io,fc,pCp fe - mo^k,ss%^k 



^S(ko)Xcf 

°P{k3)\cf"P 
jt(fcj)Ac/ , 

«P ap{kj)Xcf 



12 



™o,fc,s(c5_fe - Sp J.) + TOo,fe,p(cp_fc - slj.) + 2kkk {ss,kCs,k + sp,kCp,k) 



+ 



kkkicp^k ~ Sp,k) ~ {mo,k,S + mQ^k,p)spMCpM 
kkk{c%,k - '^'s,k) - (iTi-oM.s + mo^k,p)ss,kCs,k 

I 



's{kj)\crp 



P(kj)XcfUs 



L(kj}Xcf 



jt(fcj)Ac/ , dk])Xcf 
°P(k])Xcf"'S +as{kj)XcfOp 



(49) 



Following the procedure described in Section III.C, we 
have to solve the gap equations for the mixing angles 9s^k 
and 6'p.fc 



kkk{c%^k - ■Ss,fc) - (™o,fe,s + rno^k,p)ss,kCs,k 
kkkicp.k ~ Spk) - (mo^k.s + mo^k,p)sp^kCp^k 



= 
= 0, 



(50) 



respectively and we find immediately that these equa- 
tions imply 9s.k ~ 9p^k = ^fe- The dk values are listed 
in Table [IT] for different possible solutions and are in 
general not equal. The masses are given by mQ^k,s = 
J2n' "^ol/^Af'fcP, TOo,fc,p = Z^AfTOolajvfcP- The bare mass 
Too = Q.OOSGeV, is chosen at the intermediate value be- 
tween the mass of the up and down quarks [l0| . Finally, 
with the values kn and fc22 obtained (see Table [l| we can 
calculate the pair-energies es^kk, ep.fcfc and the ground 
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Sol 


ei 


62 




ep,ii 


£0,11 


es,22 


ep,22 


£0,22 


1 


0.76331 


0.76828 


0.18111 


0.18111 


4.34658 


0.23372 


0.23372 


5.60922 


2 


0.75749 


0.76737 


0.14343 


0.14343 


3.44223 


0.22199 


0.22199 


5.32784 


3 


0.77226 


0.77305 


0.30461 


0.30461 


7.31072 


0.32394 


0.32394 


7.77446 



TABLE II: Solutions of the BCS equation, using the a and /3 values from Table I. The lowest energy solution is comes from 
the Set 2, as can be seen by the es,ii and ep.n values. 



State energy eo,fcfc: 



12 



'kkss,kcs,k + "io,pC5 - mo^sSsj, 



£p,kk 



?'«o,s(c|_fc - Sp^fc) + rno^p{cj,,^ - j.) + 2kkk iss,kCs,k + sp,kCp,k) 

I 



£0,kk, 



(51) 



the results are shown in the Table HIl 



Taking the sohition which yields the lowest energies we 
obtain for the kinetic energy of the Hamiltonian 



+.P,n(b;^,i)6^(^*'+4(,.)d-(^i')(-^2) 

The energy of the lowest s and p orbital are degener- 
ate. This is due to the structure of the kinetic energy. A 
residual interaction in the potential term should remove 
it. Single-particle states belong to the spin color (1, 0) 
and flavor-spin i representations. The antiparticles be- 
long to the conjugate representations. In order to build 
the low lying meson states, we have to put one parti- 
cle and one anti-particle in the lowest level, with energy 
0.143 GeV. Thus the energy of the lightest mesons would 
be expected around 0.286 GeV, in good agreement with 
the scale of excitations. Considering that a baryon needs 
three quarks, in order to be in a colorless state, one ob- 
tains the first states about 0.429 GeV, which is by about 
a factor two lower than the empirical value. 

As shown in section II. A, the potential term is pro- 
portional to the second order Casimir operator of SU (3) 
color. This term commutes with the kinetic energy and, 
thus, one can diagonalize the kinetic energy indepen- 
dently to the potential term. 



V. CONCLUSIONS 

In this work we have analyzed a QCD-inspired 
Hamiltonian, considering a constant interaction between 
quarks. The constant interaction represents an average 
over all residual interactions. All orbital levels were taken 
into account in the calculations. The color part was de- 
scribed by a SU{3) group and the flavor part was de- 
scribed by a SU{2) group (only up and down quarks 
have been included in the model). Nevertheless, it can 
be easily extended to SU{3), by adding an additional 
mass-term for strange quarks to the kinetic energy. No 
gluons have been included dynamically. The system was 
confined to a finite volume, with the length of an average 
harmonic oscillator. 

The potential interaction resulted in the color-spin op- 
erator (second order Casimir operator of the color-5L/(3) 
group). This interaction divides the colorless states from 
the colored ones, shifting the colored states towards 
higher energies, being the shift proportional to Vo, the 
strength of the constant interaction. 

The kinetic energy operator could be written in terms 

of the sum of the Kj_ operators, restricting to two- and 
three-level sub-systems and excluding the mass term. We 
proved that they satisfy an SU{2) algebra, by adding 
a i^o^ operator which counts the difference between 
quarks in the upper and lower orbital levels. This sug- 
gested that an analytic or quasi-analytic solution may 
exists also for the general case. This solution was found 
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after introducing an unitary transformation, to the quark 
creation and annihilation operators, foUowed by a simple 
BCS-type of transformation. 

The main features of the calculated spectrum are: 

i) Its structure is already quite rich, and resembles the 
physical meson and baryon spectrum, 

ii) the mesons and baryons consist of two and three 
partons, respectively, which once transformed into the 
original quark operators correspond to a back-ground sea 
of many quark-antiquarks pairs. 

As shown, though QCD is highly non-pcrturbative, 
there is still a lot of room for finding analytic solutions 
for specific sectors of the theory, particularly, by applying 
algebraic and non-perturbative transformations based on 
many-body techniques. 
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Appendix A: Constructing the potential interaction: 



The potential term is given by the color-color interac- 
tion 



V = J dxdyxp\x)T-'tPix)V{\x~y\)xp\y)T'''iPiy), 

(53) 



We acknowledge financial help from DGAPA and from 
the National Research Council of Mexico (CONACyT), 



where the fermion operators xp^ and have the following 
structure 



J 



i,\x) = {xp\{x),4{x)) 



ipix) = 



iNhnacf Nlmacf 



(54) 



So the potential term becomes 



V 



^ / dxdyV{\x - y\) 

aiNi,limiji\i(TifiCiLM 

6"=('^=''^'^^^'^^^^=^M^2m2,i-72b2A2)i?JV.i.(a;)i1.m.(f^.)X.2'5/i/2'5aiai<5,i.. 



(55) 



The angular part of the double integral is given by the 
following expression: 
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I 

where the residual interaction can be written as The second expression conies from the orthogonahty 

of the Legendre polynomials and the third is a useful 
relation for the Legendre Polynomials. Therefore, the 

V{\x-y\) = J2^lPl{cos9) 

double angular integral can be separated in two single 



angular integrals: 



^Al = (^^2^) j ^d{cos6)PL{cose)V{\x-y\) 
Pdcose) 



{T^^Y.^*LM{^^)yLM{^v)- (57) 



dn.,dnyYi\,^^{n,)Yi,,„.,{n^)v{\x - v\)Yii„,,{ny)Yi,„,,{^y) 

47r 



LM 



LM 



2L+1 

/ 47r 
\2L+1 



(-1) 



nil +-A/4-m3 



C?^y^/3-m3(Oy)^/4m4(ny)^LA/(^y) 



h h L 




(-1) 



— mi 7712 —A/ 



{2h + l){2l2 + l){2L + l) 



An 



(2/3 + l)(2Z4 + l)(2i+l) 



47r 



/a I4 L 




I3 U L 
— TO3 TO4 M 



(58) 



With the above expressions , the potential term takes 
the following form 



y= E 

aiNi,limiji\i(7ifiCiLM 



{-^ )c2" XaiXa-2"flf2"aia2'>a 



ai{NiJi,^)jlXlcih^-^ 'C2 

a3(Af3,i3,^)j3A3C3/3^^ ' '^i" Xa4 O/3 /4 0a3a4 

a;pd|a;||ypd|y|i?^^,^(|a;|)i?Ar,z,(|a;|)i?^,,,3(|y|)i?w,^(|y|) 



d{cosd)PLicosd)V{\x - y\ 



2L+1 



nv(2/.+i)(2,.+i) ' ' 



(-1) 



rni+M+vi3+h~l2+h-l4+li-k-^i+h-h-^2+h-k -^3+14-^-^4 



^'^^ {liO,hO\LO) {hO,hO\LO) 



h \ ji 

mi Gi — Ai 

^3 \ J3 

ms 0-3 -A3 



i2 2 

7712 (72 — A2 

m4 (T4 — A4 



—mi m2 —M 
—771,3 "^4 M 



(59) 
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-(Tl — —(73 — ■k-\-AI — A2 — -k — A4 — i 



h \ h \ f h h 32 \(h h L 



2 2 
\ -Al mi (Tl y V -m2 A2 -(Tl I \ m2 -mi -M ) 

jl j2 i \l jl 32 L \ |.__^^i_i_^3_i_A2_i-A4-l_ 



-Al A2 -M \ I2 h h 



(60) 



Joining the resulting phase with the last three 3-J co- is obtained 
efficients of the equation Eq. , the following structure 



77130-37714 

I 33 I 



h \ 3i \ ( h h L 
m^ (T3 —A3 / V TO4 (Ta — A4 j \ —7713 m4 M 



1 

3 2 
A3 "^3 0-3 



— m4 A4 — (T3 y I 7714 — ^3 M 



\L— i-cr3—i-A2-i-A4 -l + /4 + i3+L- (3 — i—U— 7714 — 7713+0-3 



(-1) 



'3 + ^+i4+77l4+7773 — 03 



33 3i L \ I J3 J4 L 



-A3 A4 A/ ; \ /4 ^3 5 

where the 3-J coefficients can be written as 



(-l)^'^(-l) 



(61) 



31 32 L 

-Al A2 -M 

33 34 L 
-A3 A4 M 



(-1) 



jl-j2+M+ji+j2-L 



V2L+1 

I 



(.7iAi,j2-A2|L-M) 



(62) 



Another contribution to the general phase conies from 
lowering the spin and color indices of the annihilation 
operators, which gives the additional phases 



X (^_iyi-j2+M+ji+j2-L^_-^y3-j4-M+j3 +ji ~ L 

= (-l)Hj2 + i+J4(_l)M(_-^)X2-fX4^ (64) 



(^_iy2->^2+X2 (_-[^'ji4->'4-fX4 

With this the general phase is given by 



(63) 



(-1)^^ (-1)^ 



5-FA2 + 5+A4 f^_^y2-^2+X2 V4-A4+X4 



(-IF 



where (— 1)'<:'= arc the S'C/(3)-phases of color, as defined 
in i. 

Finally, the product of the SU(2) color generators can 
be written as follows: 



{Ta)ZiTlZ = ((1'0)C2, (1, l)a I (l,0)ci)i(-l)>^»((l,0)c4, (1, l)a \ (1, 0)03)1 ((1, 0) ||| T ||| (1,0))^. (65) 



The factor ((1,0) ||| T \\\ (1,0)) is a triple reduced ma- order Casimir o pcrator of S'?7(3) d which is equal to -s/S. 
trix element and its value is given by twice the second 
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The bar over an index refers to the conjugate component. 
The index 1 at the end of an SU{3) Clebsch-Gordan co- 
efficient is the multiphcity of the coupling. Using the 
symmetry properties of the SU (3) Clebsch-Gordan coef- 
ficients and changing the index a by C, we arrive at 
the expression 



3(_l)x..+xc,+xc((i^0)ci,(0.1)c2 I (l,l)C)i 

((I,0)c3,(0.1)c4 I (l,l)C)i .(66) 



Finally, the potential term is given by 



Ni,ji,Xi,h,Ci,aJ,a' J' ,L,M,C 



J 



Ql(Afl,il,|)jlAlCl/ {N2.l2,^)j2-X2-C2 



^l(A^3,.3,^b3A3C3/'^r^i,lb.-A.-c.03A3,j4-A4|iM)((l,0)c3,(0,l)c-4|(l,l^ (-1) 



(jiAi, J2 - X2\L - M)((l,0)ci, (0, 1)C2|(1, 1)C) 

M+xc 



(-1) 



22LTT / l^r«'l^llyr«'lyl^kh(l^l)^A^2/.(|a3|)i?w3;3(lyl)^A^4'4(|y| 



d{cos9)PL{cos9)V{\x - y\) 



-1 

4 



n Vi^h + l)(2j. + l)(/iO, hO\LO){kO, hO\LO) I ^1 ||{^' J 



r 



(67) 



The Clebsch-Gordan and all recouping coefficients are SU{3) Clebsch-Gordan and rccoupling coefficients, are 

calculated numerically [1, 0|- In Q a useful collection available, 
of formulas, comprising all symmetry relations known 

for these coefficients, is given. In [ll|, [T^] more recent As mentioned before, the potential term has the fol- 

versions of the computer routines, for calculating the lowing form: 

I 



V = ("i)*'(-i)'^- 

Ni,ji,Xidi,Ci,aJ,a' J' ,L,M,C 



Q(JVi,ii,|)jlAi,Ci/ (A'2,i2,5b2A2,C2 



, (ii Ai, J2A2 \LM) ((1, 0)ci, (0, 1)C2| (1, 1)C) 



/ r/ 

a J 



"'i(-'V3,'3,5)i3,A3,C3/' {N4,,l4,^)jiX4,C4 



, (jsAa, j4A4|i - M)((l, 0)c3, (0, l)c4|(l, 1)C) 



^/{2n + l){2js + l)ViN,,k,j,,L) 



(68) 



The creation operator is in a (1,0) S'/7(3)-color irrcp and 
the annihilation operator with a lower index is in a (0,1) 
S'[/(3)-color irrep. A creation operator is coupled with 
the annihilation operator to a (1,1) S'[/(3)-color irrep, de- 
noting the color generator. The quantity V{Ni,li, ji, L) 
is the intensity of each component of the interaction, and 
it reads 



viN,,k,,n,L) 



1 



2 (2L + 1) 



\x\'d\x\\y{'d\y\R*^^t^{\x\)RN2i2{\x\) 



Rkisi\y\)RN4i4i\y\) 



X / d(cos6l)Pi(cos6l)F(|a;|,|y|,cos6') 

r4 



nti Vi'^h + l)(2j. + l)(^iO, /20|£0)(;30, ?40|L0) 
V(2ji + l)(2j3 + 1) 



jl ^1 5 1 i J3 ^3 5 



(69) 
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Taking into account that the only allowed value of L 
is L=0 and using a constant potential Vq as residual in- 
teraction, implies the following simplifications: 



d{cos9)PL{cos9)V{\x\, \y\,cose) = 2Vo 
x\^x\\yfd\y\R*^^i^i\x\)RN,i,{\x\) 

Rkhi\y\)RN.u{\y\) = 



\xfd\x\R%^,^{\x\)RN,i,i\x\) 
y\'d\y\R*^^,MRN,u{\y\) = 



(/iO,Z20|LO) = 



V2/i + 1 



Ji h 

^1 '2 5 J v/(2.7i + l)(2^i + l) 

h U \ ] ^(2j3 + l) (2^3 + 1) 
therefore the intensities ^(A^i, ^j, Ji, -Z^) are given by 



Appendix B: Constructing the kinetic energy 
operator: 

We start from the expression of the kinetic energy term 
obtained in section II, without the mass term. The con- 
tribution of the mass term will be considered Appendix 
C. 



00 N+l 

JV=j + i JV' = W-1 Ac/ 
00 N+l 

7V=j>i Af'=Ar-l Ac/ 



(71) 



where N varies in steps of 2, and the kjNN' constant is 
given by 



i J r^dr 



R 



jv'0-i)(0 (72) 



We use the 3-dimensional harmonic oscillator, i.e., 
{Rn'i{r) = iV„'/exp(-^)r'L^t^(7r2)) as radial func- 
tions and the following convention: N' = ^^-j-^, which 
is the conventional notation for the harmonic oscillator 
(at the end we will come back to the TVj-notation, where 
I ~ j — and I' = J + ^)- The factors kjjsiN' can be 
written now as 



Klnn' 



r dr 



(73) 



The result of applying of the operator — on the 
radial function is 



-7r'+iexp(-^)L'+^(7r2)} . (74) 

From the recurrence relations of the Lagucrre polyno- 
mials 



d 

dx 



(75) 



as 



the first relation the equation Eq. (|74p can be written 

Nr,i{r'eM--^)i-rL'n'hir^) 
-7r'+iexp(-2I^)(<t^+*(7r2) - Llt'^XHir^))} 



(76) 



and 
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r'exp(-2l!)(-27r)<t!t"(^) ■ 
We obtain the final expression 



(77) 



n = n - l^N' = N + l 
n = n' => N' = N -1. 



(80) 



N, 



.'+iexp(-2l!)<tl+*(7r2) 



n'l I -27? 

-7r'+icxp(-^)<r+^(7r^) 
V+^exp(-2l!)Litlt*(7^')} 



= Nn'i -7 



(78) 



Therefore, the connections for the harmonic osciUator 
are given by 



The quotient to the normahzation constants is given 



by 



Nn',l , , . Nn',l ^ ^ 
-11— On,n' -lOv ,1+1 - IJ— d„,„'d/'J+i. 



N, 



n',l+l 



-Sn,n'-lSl' .1 + 1 



AT — L'^l' ,1 + 1 



2(»')! 



r(ri'+i+j 



2 111 
'74^2 



r(K'-i+;+i+f ) 



2(n'-l)! 



3 I ' + 1 
/ 4 "t" 2 



2^^^ 0N',N+1 



(81) 



and 



(79) 



2(n' 



r(n'+i+|) 



2 3 I i 

0/4^2 



r(n'+i+l+|) 



2(n')! 



In order to translate this expressions to the original 
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■5n'.n-i 



(82) 



N'-l 



and 



so for the 



notation we write n' — — 2 — "^^^ — — 9 — i 

first and second term of Eq. [79] one has the following The kinetic energy term is written as K = K^^^^^ 
selection rules (3+4) ■^yii;]-^ 



_^(l+2) 



E E 



Af+l 

E E 

j N=j+^ N'=N-1 \cf 



^"i(JVj + i,i)jAc/" 




+ 6^ 



i(Wj + i4)jAc/ 



b2 



i(W'j-i.i)jAc/ 



(83) 



00 Af+l 

-EE EE 

J Ar=j+i 7V'=Ar-l Ac/ 



'iV- 



V71 



-z ON'.N+i + V7\/ dN',N~-i 



r^(^'j-i^)jAc/'' +b_i(jv',,-i, lb-Ac/'' j 



(84) 



Finally, the operators K+ = A'^^^ + A'^'^' and K- = 
^(2) _|_ ^(4) i^g^yQ structure 
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K. 



oo N-1 
3 N=j+^ N'=N+1 \cf 



On',N+i + 



t: On'm-i 



oo N+1 

HV7)E E E E 

J N=j+i N' = N-1 Xcf 



N',N-1 



i(Arj + i,ibAc/ 



(85) 



Now WG redefine, for each set N,j with I = J + ^, 
the fermion creation operators multiplying them by {—i) 
and the corresponding annihilation operator by In 
this way, the anti-commutation relations are preserved. 
Explicitly, the mapping is given by 



(86) 



±i(Afj + i,i)jAc/ 



"±l(NJ+i,i)j\cf 



With this the kinetic energy parts acquire the form 



oo N+1 

(V7)E E E E 

] N=j+^ N'=N-1 Xcf 



'n-j + I n + j + I 

:z ON', N+1 + \ t; dN',N-i 



oo N+1 



K = 



(x/7)E E E E 

'J N=j+^ N'=N-1 Xcf 



dN\N+l + \ 7^ dN\N-l 



Appendix C: Contribution of the mass term 

Substituting the fermion fields into the expression of 
the mass term, we obtain 



(87) 



J 



dx'il)\x)j3moip{x) ^ J dx yi{)\{x) , xp^i^ 



-mol J \ tp2{x) 



^ 2 1 

dx ^ -Cri|jlAl)(;3m3, -CTgljgAg) 

NiNshlsmimsjijaXiXsaicracf 

R*^^,^i\x\)RNaU\x\)V^^,{n,)Yi,^A^,^^^^^ 

^ 2 1 

dx ^ {hm2,-(J2\j2>^2){km4,-^0-4\j4>^4:) 

N2Nil2lim2mij2jiX2Xia2 0-40/ 



21 



I 

The radial and the angular integral leads to the restric- plying also the mapping of the fcrmion creation and an- 
tion SNiNaSi^i^SmimsSaiaa and (JTVsATi^ijii^msmi^cracri , re- nihilation operators, as given in Eq. ^^): 
spectively. Thus, the mass term can be written as (ap- 



Nl]\cf 



] Xcf N=3 + ^ 



where N and N' vary in steps of 2. Now, we apply the 
same unitary transformation as described in the text, i.e.. 



"±i(Ar'j-i,i)iAc/ 



Q^iVfc^iji j-(fc 7+j- -)j\cf With this, the mass term is written 



dxxjj\x)f3moilj{x) ~ 



EEE -0, 

j Ac/ q ^ 



EEE 

j Xcf k 



fe-i(fcj + i4)jAc/^ 



(91) 



Due to the method described in section III the only re- 
maining terms correspond to k=q. With this the masses 
become 



(92) 



n-l 



Finally, the mass contribution to the BCS kinetic energy 
is given by 



K 



BCS 



EE{K^-., 

Ac/ kj 



7,t /,(fej)Ac/ 
^-i(fej)Ac/"j-i 



jt(fej)Ac/ , 
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IS- 1 



12 



Lt jt(fci)Ac/ , ,{kj)Xcf 

Vi(fcj)Ac/^-i +rfj-i(fej)Ac/fj+i 



APPENDIX D: The kinetic energy term for three 
levels and a fixed and j 



(93) 



The Dirac picture of the three-energy levels implies six 
states, three with positive energy and three with negative 
energy. Fixing j and N, we are working in a determined 
column with total spin j and the possible connections 
{N - 1, J - i) ^ {N,j + i) and {N,j + I) ^ {N + 
liJ ^ 5)- The operators and K of this case have 
the following structure, respectively, 



^(JV,i+i,i)jAc/" 



"i(W+lj-i,ibAc/ 



"i(JV-lj-i,i)jAc/ 



N+j + i 



N-3 + i 



-|(Wj + i,ibAc/" 



(W+lj-i,ibAc/ 



55Wi+ii)jAc/ 



5i(Afj + i,i)jAc/ 



(94) 



The commutation relation of these operators is given 
by (omitting the sub-indices jXcf for simplicity) 



+ 





2 




N 


-j^ 


3 




2 




N 


-j^ 


_ 3 




2 




N 




3 




2 




N 




3 




2 




N 


- 


3 




2 




N 


- j4 


_ 3 




2 




N 


+ j4 


3 



AT 



N + J + 



N + j + l 



"-1(^+1,^-1,1)" 



AT 



^-i(iV,, + i,i) 



"1(^-1,^-1,1)" 



TV 



i(7V-l,j-i,i) 



N 



23 
(95) 



[ Ko ,K+] = 
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N-j + l 



N-j + l 



+ 2 



N + J + 



^[2iV + 3]ii:+ 



(96) 



Defining 



Choosing 



K± 

we arrive at the relation 



= VK± 



(97) 



(98) 



7 



(2A^ + 3) 



(99) 



we arrive at the weh known SU{2) algebra. Of course, 
this has to be repeated for whieh is easily obtained 
by taking the adjoint of K^. 
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